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Dynamics of Skyrmionic spin texture in the spin-1 Bose-Einstein condensate (BEC) is examined by 
analytical and numerical means. We show the Skyrmion (coreless vortex) to be inherently unstable 
in the sense that the state initially prepared purely within the anti-ferromagnetic (ferromagnetic) 
order parameter manifold inevitably evolves into a mixture of both. The vorticity-dependent drift 
in the presence of the trapping potential also contributes to the disintegration of the initial spin 
texture manifold. We argue that the notion of a Skyrmion as a topologically protected entity 
becomes ill-defined during the dynamical evolution process. 
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I. INTRODUCTION 

Optical traps enable the possibility to host more than 
one hyperfine spin state of cold atoms simultaneously, 
therefore opening up the popular study of spinor Bose- 
Einstein condensate (BEC) since the first experimental 
realization [1]. The pioneering theories of the spin-1 BEC 
worked out by Ho [2], and independently by Ohmi and 
Machida [3], allowed the understanding of the ground 
state structures as either anti-ferromagnetic (AFM) or 
ferromagnetic (FM), and elementary excitations. Besides 
the equilibrium properties, the dynamical behaviors also 
received great experimental attention [4] later on, breed- 
ing new research interests. 

The concept of Skyrmion comes from a model for 
baryons in nuclear physics [5], and has been regarded as 
a topological particle of great significance in condensed 
matter physics as well. In AFM spinor BECs we can 
also predict the existence of the counterpart which cor- 
responds to a metastable excitation, while in FM spinor 
condensate the phrase "coreless vortex" is adopted some- 
times [6, 7]. Bogoliubov-de Gennes equations are very 
useful to study the stability of such spin textured states in 
cold atom systems [8, 9]. Thermal fluctuations, etc, can 
limit the Skyrmion lifetime which is the key element to 
make sure we can observe, even manipulate Skyrmions in 
BEC system. Previous experiments have already success- 
fully created the Skyrmion in both AFM (^^Na) [10, 11] 
and FM (^^Rb) [12] condensates. The decay process into 
the ground state was also observed [11]. However, few 
discussions have been made on the exact temporal dy- 
namics of Skyrmions before reaching the thermal equi- 
librium. Here we focus on the pure spin-1 BEC system 
at zero temperature, and study the time evolution behav- 
iors starting from either AFM or FM states supporting 
the Skyrmionic spin texture. 

The paper is organized as follows. In Sec. II, we intro- 



duce the concept of Skyrmion (coreless vortex) in the 
spin-1 condensate for both AFM and FM cases, and 
in Sec. HI we analytically discuss the possible decay of 
Skyrmion (coreless vortex) originated from the break- 
down of AFM (FM) phase in the dynamical process. Nu- 
merical simulations of the real time evolution behaviors 
are made in Sec. IV to show the detailed decay phenom- 
ena. Three different initial configurations are considered, 
exhibiting different evolution patterns. Finally in Sec. V, 
we summary our results and make corresponding discus- 
sions. 



II. GROUND STATE AND SKYRMION 

In the mean-field approximation, the energy functional 
for the spin-1 condensate without the trapping potential 
is commonly given by [2, 3] 

= £ y"(V*t).(v*) + i Jp^[co + c,{Sn{l) 

where cq = (ao + 2a2)47rfi^/m and C2 = (^2 — ao)47rfi^/m 
characterize the interaction strengths, {i = 0, 2) are 
the s-wave scattering length in the two-atom i-th scat- 
tering channel, and m is the particle mass. The spin-1 
condensate wave function can be decomposed as 
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where p and are the total density and overall phase, 
and S = rj^Frj describes the spin vector with rj the unit- 
modulus spinor field obeying rj^rj = 1. The 3x3 angular 
momentum operator for spin-1 condensate is represented 
hy F = {Fx, Fy^Fz}. It is a common practice to dis- 
tinguish the two regimes of the ground state according 
to C2 being repulsive or attractive. The former corre- 
sponds to the AFM case where the spin average becomes 
zero, |S| = 0, to minimize the spin-dependent interaction 
energy, whereas the latter is the FM case and the spin 
average is maximized, |S| = 1. With the aid of spin ro- 



tation operator U{a, /3, 7) 
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a, /3, and 7 are the Euler angles, we may express the two 
fields as [2, 3] 




(3) 



Throughout the paper we denote rji and 772 for AFM 
and FM order parameter manifolds, respectively. The 
two angular variables a and f3 can be used to con- 
struct a unit vector d = (sin /3 cos a, sin /3 sin a, cos /3) 
which in turn can give rise to a topological spin tex- 
ture of the Skyrmion (coreless vortex), where the integer 
Q = J dxdyd • {dxd x dyd)/{4:7r) denotes the topologi- 
cal charge. Typically, Skyrmions refer to a configuration 
where the Euler angle a corresponds to the azimuthal 
angle (j) of the plane, and /3 is a function of radial dis- 
tance r varying from at the origin, to tt at infinity or 
the boundary of the condensate. Skyrmions can be im- 
printed in spinor BECs by a spin rotation method. In a 
recent experiment [11], the Skyrmion created in the AFM 
^^Na condensate decays into the ground state over time, 
producing what appears to be a half-quantum vortex- 
anti- vortex pair [13] within the condensate in the inter- 
mediate phase. Motivated by the dichotomy of the (the- 
oretically) expected topological stability of the Skyrmion 
in the AFM BEC and its smooth decay found in the 
experiment, it is a timely exercise to carry out a more 
critical analysis of the Skyrmion dynamics in the spin-1 
condensate. 



III. EQUATION OF MOTION ANALYSIS 

Some "pathology" in the spin-1 AFM dynamics was in 
fact noted early on [14]. The kinetic energy obtained for 
the AFM wave function ^1 reads 



El 
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while the action part Si = iJ^jSt^i becomes Si = 
— / pdfO. As it happens, there is no term in the action 
responsible for the dynamics of the d- vector! By contrast 
the FM spinor wave function ^2 = v^^*^^2 gives rise to 
S2 = — J pdtO + J (cos /3 — l)dta, where the second term. 
Berry phase action, is responsible for the spin dynamics 
of d governed by Landau-Lifshitz equation of motion. 
The discrepancy is further illustrated by the examination 
of the overlap integral between adjacently located wave 
functions ^(r) and ^(r — 6t) for the two cases [15, 16], 



<^r)«'i(r) 
<5r)*2(r) 



p(r)e^ cos/3(VQ:-5r)^ 



The overlap of the adjacent FM wave functions produces 
the Berry phase factor, which is absent for AFM wave 
function overlap. 

A way to cure the pathology of the d- vector dynam- 
ics in the AFM manifold was suggested by the authors 
of Ref. [14], who considered small fluctuations away 
from the AFM manifold and obtained, by integrating 
out the fluctuations, an eff"ective action for d that is 
quadratic in time ~ (9d/9t)^, rather than flrst-order 
as expected in superfluid vortices [15, 16] and magnetic 
Skyrmions [17, 18]. A similar idea was explored by 
Ruostekoski and Anglin, who numerically observed the 
spontaneous deformation of the monopole core into an 
extended defect called the Alice string, and attributed 
the phenomenon to the energetic balance of AFM and 
non-AFM components in the wave function [19]. 

Here we want to offer another perspective of this dy- 
namical problem. When we deflne the spinor flelds in 
Eq. (3), we rotated the two bases (0, 1, 0)^ and (1, 0, 0)^, 
but ignored the third one (0, 0, 1)^. Euler rotation of this 
third basis yields 



(6) 



which also corresponds to the FM manifold. Together 
the three flelds rjj (j = 1, 2, 3) form a complete, orthogo- 
nal set obeying rj^jTjk = Sjk- An arbitrary spinor rj must 
therefore be a linear combination of the three basis vec- 
tors. Now we consider the initial AFM or FM state, and 
let it evolve for a small time step At to see the compo- 
sition of the intermediate state, i.e., ^(At) — ^(0) ^ 
-iAtH{0)^{0), where H = -h^V'^/{2m) + cq^^^ + 
C2(^tF^) -F. 

Firstly, for the AFM state we set ^(0) = ryi, and the 
action of Hamiltonian on it yields 
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with 



«i = + (diaf sin^ 13] + cq, 
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{Ai+A2 + A3-A4)e'\ 



'2m 

-^{-Ai+A2 + A3 + A4)e-'\ 
2m 
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where we have deflned Ai = {dia^ sin P cos /3/a/2, ^2 = 
i{^fa)smp/^f2, = i^/2{dia){diP) cos and A4 = 
df(3/V2. Note that the Einstein convention is used, and 
i = x^y for two dimensions. In general, Uj {j = 1,2,3) 
are non-zero. As a result, the initial pure AFM state 
would evolve into the mixture containing the FM com- 
ponent. The disappearance of C2 in uj indicates that C2 
plays no role in the beginning of the evolution. 
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For the FM initial state, ^(0) = r/2, we have 



Hr]2 = ~^^^^2 + Co 7/2 + C2d • Fr]2 = Vjr]j, (9) 



with 

+ ^ {dia)'^ sin (3 cos fSje''^ , 
v2 

^2 = - ^2 + ^3 - ^4 - ^5 - ^e) 

2m 

+Co + C2, 

^3 = + ^2 + ^3 + ^4 - ^5 + ^6)6-"^^ 

2m 

(10) 

where we have Bi = [(Sia)^ + (dij)^ + ^9?Ci^]/2, ^2 = 
[2(S,a)(a,7) + iah](cos/3)/2, ^3 = sin^, 
54 = + i5f7 + (9.7)']/2, ^5 = [2(S,c.)(S,7) + 

idfa]{cosf3)/2, and Be = (9^a)^(cos^ /3)/2. Immediately, 
we can draw the conclusion that the initial pure FM state 
would also be mixed with AFM component during the 
time evolution. Additionally, positive (negative) C2 may 
enhance (reduce) the effects of cq, which will be further 
confirmed later. 

The main lesson of the above short-time analysis is that 
we can readily deduce the source of AFM-FM mixing to 
be mainly the inhomogeneity in a and /3, that is to say, 
the kinetic energy carried by the inhomogeneous initial 
condensate wave function. It is the spin texture in the 
initial configuration which drives the mixing. A textured 
localized object such as a Skyrmion (coreless vortex) in 
a pure AFM (FM) condensate inevitably evolves into a 
mixed phase due to its intrinsic inhomogeneity. 



IV. DESCRIPTION OF NUMERICS 



general consideration we introduce the displacement of 
the density peak from the trap center by ^p, while 
is the displacement of Skyrmionic spin texture. We as- 
sume both and are small quantities (compared to 
the condensate size) in the x direction, i.e., = (^pe^^, 
and = (,s^x' The density distribution takes the Gaus- 
sian profile p{r — ^p) = e"^'^"^^^ /\/27r which is the single 
particle ground state in the harmonic trapping potential. 
Our choice of the density profile is not the exact many- 
body ground (or metastable) state of the given Hamil- 
tonian. For the sake of evolution, derivation from such 
exact state is desired, therefore we choose the current 
form for simplicity. For the AFM spinor fields we have 



cos/3(r - Is) 
V 75e^^^''"^^^sin/3(r-|s) 



(12) 



where (/)(r — |s) is the azimuthal angle, and cos /3(r — |s) = 
[(r-^,)2-l]/[(r-C,)2 + l],sin/3(r-C3) = 2|r-C3|/[(r- 
|s)^ + 1] is used. For the FM case, in view of the recent 
experiment [12], we may adopt 
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_i^e^^('-^s)sin^(r-|s) 
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where /3(r — |s) = 7r(l — e~l^~^^l) varies from to tt. 
Note that similar results are always obtained regardless 
of the choice of the exact radial dependence of /3. We may 
define ^ = (t/^+i, '0o, V^-i)^, where ipi (i = +1, 0, —1) de- 
note the wavefunctions in the corresponding hyperfine 
states. Initially, the winding numbers of '^o, and 
for AFM and FM cases are (-1,0,1) and (0,1,2), 
respectively. In the following calculations, the real time 
evolution method is used, and we set cq = 50.0. Con- 
sidering the possible experimental setup, we discuss the 
following three kinds of situations with regard to different 
choices of and ^s- 



Exact time evolution behavior can be obtained by nu- 
merically solving the full set of Gross-Pitaevskii (GP) 
equations in the presence of a trap 



'9 ^ 



2m 



V + co(^^^) + C2(^^F^) • ^, 



(11) 

where V = tykjo'^y'^ /2 (r = xe^ + y^y) is the 2D har- 
monic trapping potential. To make all quantities di- 
mensionle ss we s et the energy, length, and time scales 
as /kj, ^Jh/muj^ and l/o;, respectively. As we are 
mainly interested in the Skyrmion (coreless vortex) dy- 
namics in the spin-1 condensate, initial st ates are cho- 
sen to be eithe r ^AFM(r ,0) = ^^ p{v - ip)'qi{v - i^) or 
^FM(r, 0) = \/ p{v - ^p)ri2{r - |s) which represent the 
Skyrmion (coreless vortex) in pure AFM (FM) state. For 



A. = ^3 = 

In this case, both the density and the Skyrmionic tex- 
ture are centered, enabling the rotational symmetry to 
be preserved during the whole evolution process. As dis- 
cussed in Sec. Ill, either a pure AFM or a pure FM initial 
state will evolve into a mixture of both components. In- 
deed, as shown in Fig. 1(a), non-zero |S(r, t)| {t > 0) 
indicates the appearance of the FM feature out of the 
AFM initial state, while in Fig. 1(b), the deviation of 
|S(r, t)| from 1 shows the loss of FM character although 
the initial state was fully ferromagnetic. 

Naturally it becomes problematic to give the precise 
meaning to the Skyrmion now that the very manifold 
(either AFM or FM) on which it was defined gets disinte- 
grated over time. It might be possible, through decompo- 
sition of the condensate wave function in the basis space 
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spanned by r]j (j = 1, 2, 3) as introduced in the previous 
section, to deduce the d(r, t)-vector and the Skyrmion 
density d • {dxd x 9yd)/47r formally at any given time t 
for any wave function ^(r, t). This is a subtle and impor- 
tant issue which we will treat in a later publication. In 
the mixed AFM and FM state, the winding numbers of 
i^i = +I7 — 1) are still preserved during the evolution, 
and the vortex (anti-vortex) cores always coincide with 
each other. Since some Skyrmionic features continue to 
be maintained over time, we may name this structure in 
the mixed state as a "pseduo- Skyrmion" for the conve- 
nience of description. 



|S| 0.5 



(a) AFM 



(b) FM 





FIG. 1. Radial dependence of \S{x,y = 0,t)| for the initial 
Skyrmionic configuration defined on (a) AFM and (b) FM 
manifolds taken at time t = 2.0 with C2 = 0.0. Other values 
of C2 from -20 to 20 were tried, supporting similar results. 



B. = and > 

In the previous subsection the initial Skyrmionic con- 
figuration obeyed the rotational symmetry. The sub- 
sequent disintegration of the AFM or FM component 
was largely due to the inherent dynamics of the con- 
densate, unrelated to effects arising from the trap po- 
tential. Now we consider the case when the Skyrmion 
(coreless vortex) is initially displaced from the trap cen- 
ter. In the spin rotation method as adopted in Ref. [11], 
by tuning the position of the zero-field center of the 3D 
quadrupole magnetic field, we may manipulate the po- 
sition of the Skyrmion (coreless vortex). Therefore, our 
setup is experimentally feasible. As a comparison, for a 
single component condensate it is well known that a dis- 
placed vortex would precess around the trap core driven 
by a radial buoyant force and a gyroscopic Magnus force 
together [20]. 

For our spin-1 AFM initial state as shown in Fig. 2(a), 
the anti- vortex of and vortex of ip-i are located 
at the same position ^s- However, they feel the oppo- 
site forces, and move in opposite directions as shown in 
Fig. 3(a)-(f). Note that similar behavior in the displaced 
AFM monopole dynamics of spin-1 condensate was dis- 
cussed in Ref. [21]. The splitting of the vortex and the 
anti- vortex cores implies that description of the conden- 
sate wave function based purely on AFM manifold is no 
longer tenable. To quantify the breakdown of the AFM 
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FIG. 2. (color online) Initial density and phase profiles for a 
displaced Skyrmion within the AFM manifold (a), and those 
of a displaced coreless vortex within the FM manifold (b). 
The first, second, and third rows correspond to V^o, and 
'0-1, respectively. Initial displacements are both at = 0.5. 
The plus signs (+) mark the positions of vortex (anti-vortex) 
cores. 



phase, we define 



/afm(^) 



/dxd^||V;+i(r,t)p-|V;-i(r,t)|- 
/dxd^||V^+i(r,t)P + |V._i(r,t)|^ 



(14) 



as a measure of the loss of AFM component. Pure AFM 
condensate would obey |?/^+i(r, t)p — |?/^_i(r, t)p| = 0. 
In Fig. 3(g) the obvious trend of increase of /afm(^) 
from zero confirms our picture that vortex- anti- vortex 
core separation implies the loss of AFM condensate. As 
t ^ 0, the effects of C2 caused by its value and sign are 
reduced as clear from the figure, confirming the analysis 
in Sec. III. 

The simulation results suggest another way for the ini- 
tial Skyrmion structure to disintegrate. The inherent 
dynamics of the vortex in the superfluid follows that of 
electrons in quantized magnetic field [15, 16]. The trap- 
ping potential provides the confining force, the analogue 
of electric field, for the vortex so that the drift motion 
orthogonal to the confining force occurs. The charge of 
the vortex is opposite to that of the anti-vortex, so the 
drifts occur in opposite directions. Unless there exists an 
enormous force binding the vortex and the anti-vortex 
cores together, the drift mechanism will inevitably sep- 
arate them in opposite directions and results in disinte- 
gration of the Skyrmion. 

As for the FM initial state in Eq. (13), both vortices 
in tjjo and ip-i have the same sign of the vorticity as seen 
from Fig. 2(b). Therefore, their cores move in the same 
direction when initially displaced from the trap center, as 
shown in Fig. 4(a)-(f). During the evolution we observe 
that the / = 2 vortex in t/j-i is unstable and decays into 
two / = 1 vortices, as shown in Fig. 4(c) and 4(f) [8, 9]. 
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FIG. 3. (color online) Real time evolution starting from a 
displaced Skyrmion within the AFM initial manifold as shown 
in Fig. 2(a) for C2 = 0.0. (a)-(c) Snapshots of densities of 
'0+1, '00, and ip-i, respectively, at t = 2.0. (d)-(f) are the 
corresponding phases. The plus signs (+) mark the positions 
of vortex (anti- vortex) cores, (g) Time dependence of /AFM(t) 
for different values of C2. 
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FIG. 4. (color online) Real time evolution starting from a 
displaced coreless vortex within the FM initial manifold as 
shown in Fig. 2(b) for C2 = 0.0. (a)-(c) Snapshots of densities 
of '0+1, '00, '0-1, respectively, at t = 2.0. (d)-(f) are the 
corresponding phases. The plus signs (+) mark the positions 
of vortex cores, (g) Time dependence of /FM(t) for different 
values of C2. 



To characterize the decay of the FM feature we define [22] 



fFuit)- 



Jdxdy{\i,+i{r,t)\ + \i,.,{r,t)\-^f 
' J dxdy{\i,+,{r,t)\ + \i>-i{v,t)\ + 



.(15) 



The definition is motivated by the observation that, in 
perfect FM state, |'0+i(r, t)| + |'0_i(r, t)| = is guar- 
anteed. From Fig. 4(g) we can observe the obvious trend 
of increase of /fm(^) from zero. Another conclusion from 
the figure is that the positive C2 may enhance the stiffness 
of the initial FM manifold, while the negative C2, on the 
contrary, would favor faster speed of the disintegration 
from the FM state. 



C. = 6 > 

Now we consider the case when the center of the den- 
sity is displaced together with that of the Skyrmionic 
texture. This situation is also realizable in experiments. 
Firstly, we create a Skyrmion (coreless vortex) centered 
in the harmonic trap following the standard procedure. 
After that we suddenly shift the trapping potential, 
therefore both the density and the Skyrmion (coreless 
vortex) core would be displaced with respect to the new 
trap center. During the time evolution the whole con- 
densate moves like a pendulum about the trap center, 
accompanied by the disintegration of initial AFM or FM 
states. Contrary to the splitting of vortex- anti- vortex 
pair or the breakdown of / = 2 vortex seen in the previous 
subsection, now the pseudo- Skyrmion structure is always 
protected throughout the evolution for both two initial 
configurations. We show in Fig. 5 the time dependence of 
the displacement of the pseudo- Skyrmion from the trap 



center, where the periodic behavior with only one pe- 
riod of oscillation is presented. We may approximately 
describe the periodic motion of the pseudo- Skyrmion as 
the semi-classical Newtonian equation of motion [14]. We 
conclude that here the dynamics of the density p domi- 
nates over that of the spin texture and effectively acts as 
a binding force tying the vortex and anti- vortex cores to- 
gether. Due to the nonlinear interactions, deviation from 
the perfect periodicity is observed. 




FIG. 5. Time dependence of the pseudo-Skyrmion displace- 
ment from the trap center for both AFM and FM initial 
states. The difference between the two is too small to dis- 
tinguish here. Initial displacements are = = 0.5, and 
C2 = 0.0. Other values of ci from -20 to 20 lead to similar 
quasi-periodic behaviors. 
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V. SUMMARY AND DISCUSSION 



Motivated by the recent experiment of Skyrmion cre- 
ation and its evolution in the AFM spin-1 condensate of 
^^Na [11], we have carried out a detailed study of the 
Skyrmion (coreless vortex) dynamics based on the GP 
equation analysis. Contrary to the standard notion of a 
Skyrmion as a metastable topological object in, e.g. non- 
linear cr- model [23], it is inherently unstable in the BEC 
condensate due to the dynamical mixing of the AFM 
and FM components during the time evolution. The in- 
stability has a different character from that expected in 
thermal or quantum tunneling of the topological object 
out of the metastable minimum ("fate of the false vac- 
uum" ) [24] . Interaction effects are also examined during 
the early stage of time evolution, exhibiting the explicit 
role of C2 in the FM case and the independence of C2 in 
the AFM counterpart. 

Numerically, three feasible initial configurations have 
been investigated for their dynamical evolution behav- 
iors, (i) For the Skyrmion (coreless vortex) initially pre- 
pared in a pure AFM (FM) manifold and located at the 
trap center, the mixing with the FM (AFM) component 
is observed. The rotational symmetric structure is pre- 
served over time, (ii) When we shift the Skyrmion while 



keeping the density peaked at the center, splitting of the 
vortex and anti- vortex centers is observed for AFM ini- 
tial state. Note that the AFM manifold breaks down 
at t > 0. For the FM initial coreless vortex state, the 
/ = 2 vortex breaks up into two / = 1 vortices, while 
rotating around the trap center in the same direction. 
This also contributes to the breakdown of the FM man- 
ifold, (iii) If both the density peak and the Skyrmion 
(coreless vortex) are displaced by the same amount, the 
whole condensate including the pseudo- Skyrmion oscil- 
lates through the trap center in analogy to the pendulum 
motion. The dynamics is governed by that of the density 
oscillation as a whole, not by those of individual vortices. 

Civen that our results rely on the energy-conserving 
simulation, they could not be applied directly to explain 
the transition into the ground state in Ref. [11]. How- 
ever, before the system reaches its thermal equilibrium, 
short time behaviors as discussed in our paper may be 
observable in experiments. 
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